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Abstract 

We present a new method for computing the wave function in the presence of constraints. As an 
explicit example we compute the wave function for the many electrons problem in coupled metallic 
rings in the presence of external magnetic fluxes. For equal fluxes and an even number of electrons 
the constraints enforce a wave function with a vanishing total momentum and a large persistent 
current and magnetization in contrast to the odd number of electrons where at finite temperatures 
the current is suppressed. We propose that the even-odd property can be verified by measuring 
the magnetization as a function of a varying gate voltage coupled to the rings. By reversing the 
flux in one of the ring the current and magnetization vanish in both rings; this can be used as a 
non-local control device. 
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The electronic wave function in quantum nanosystems at low temperatures is sensitive 
to interactions and topology such as the genus number g [1,2] (the number of holes on a 
closed surface). As a result, the wave function has to satisfy certain constraints, which 
generate conserved currents [3,4]. The implementation of the constraints is a non-trivial 
task in Quantum Mechanics [4] . The root of the difficulty is that for a given constraint the 
hermitian conjugate constraint operator might not be a constraint, therefore a reduction of 
the phase space is not possible [4] . This problem is solved by including non-physical ghost 
fields [4]. In Classical Mechanics second class constraints [4] are solved by replacing the 
Poisson brackets by the Dirac bracket and quantization is performed according to the Dirac 
correspondence principle [4,5] with the unpleasant feature that the quantum representation 
for the operators might not always be possible. Here we will solve the constraints without 
the need to introduce non-physical operators. 

As a model problem we will consider the Aharonov-Bohm geometry [6-11] for the case 
that the genus is g — 2. This corresponds to a double ring structure perfectly glued at 
one point to form a character "8" structure (see Fig. 1). Such a structure gives rise to an 
interesting Quantum Mechanical problem [12]. Gluing the two rings at the common point 
x = gives rise to a constraint problem, which was solved numerically using the Dirac 
brackets [4,12]. 

In Sec. II we present the proposed method for computing the wave function with 
constraints. We will work with a folded geometry, therefore the problem will be equiv- 
alent to a two-component spinor on a single ring. The constraint is such that at the 
common point x = 0, L the annihilation operators are identified as a single operator 
[C 1 (x)-C 2 (x)]\ x=0tL = 0. 

Using the Dirac method [3] we compute the Noether currents [4] which allow us to iden- 
tify the constraint currents. In the presence of external fluxes the constraints are translated 
into a set of equations for the wave function. The constraints induce correlations between 
the different components of the wave function. For non-interacting electrons the wave func- 
tion for N electrons is given by the Slater determinant of the single particle states, but 
the current is the same if we sum over the single electronic states. For the present problem 
we must work with the many-body wave function of the two rings (which is not a simple 
product of the two rings). 
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In Sec. Ill we discuss the constraint method with the scattering theory and show that the 
strong coupling limit between the rings corresponds to the constraint problem considered in 
Sec. I. In Sec. IV we use the constraint method to compute the many-body wave function 
for two rings in the presence of constraints. We find that the many particle wave function 
built from the single particle wave function which obeys the constraints is different from 
the many-body wave function which obeys the constraints. We show that the constraints 
impose additional relations between the amplitudes of the many-body wave function. In 
Sec. V we present the modification needed in order to include the physical geometry of the 
rings, e.g. finite thickness. As a concrete example we choose two narrow cylinders which are 
in contact on the line x — 0. Section VI is devoted to discussions. 



I. THE CONSTRAINT METHOD FOR TWO RINGS 



We consider two rings threaded by a magnetic flux $ a , where a = 1,2 represent the 
index for each ring ip a = 2tt(^^) = 2%^- = 27i0 a . The rings have a common point at 
y = (Fig. 1). The first ring is restricted to the region < y < L with the single particle 
creation and annihilation operator obeying periodic boundary conditions C(y + L) = C(y) 
and C^(y + L) — C\y). The second ring is restricted to — L < y < with similar boundary 
conditions C(y — L) = C(y) and C\y — L) = C^(y). We introduce two sets of operators: 
for the first ring < y < L, Ci(x) — C(x) = C(y), and C\(x) = C^(x) = C^(y) and for the 
second ring -L < y < 0, C 2 (x) = C(-x) = C(y), C\(x) = C*(-x) = C\y). Due to the 
folding, two equal fluxes (pi — ip 2 = (f will be described by two opposite fluxes. 

H = t dx ^°' {x) {~ Wx " T 01 ) 2 ° l{x) + ^L C * {x) { idx ~ T^) 2 ° 2{x 



i) 



The Continuity Constraint 

The Hamiltonian given in Eq. 1 will be investigated under the condition that the annihi- 
lation operators at the contact point must be identified as one operator. This is implemented 
with the help of the constraint operator: 

n = [C^x) - C 2 (x)]\ x=QtL ; V \x, N >= (2) 
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Following Dirac [3], the constraints and the time derivative of the constraints must be 
satisfied at any time. We must have r)\x, N >= and N >= 0. 

The Eigenvalue Constraint 

In order to satisfy the constraint at any time we need to show that Jr^lx, N >= 0. This 
new equation is determined by the time evolution of the constraint operator n. Following 
Dirac [3] we introduce a Lagrange multiplier A and replace the Hamiltonian H by the total 
Hamiltonian [3]: H t = H + \rfrj, where n is the constraint and A is the Lagrange multiplier. 
Using the Heisenberg equation of motion we obtain: 

i V \ X , N >= H t ]\ X , N >= ±([n, H] + \ v )\ x , N >= H]\ X , N > 0. 
In obtaining this result we have used the relations : [rj, \rfrj\\x, N >= Xr]\ X , N >. Using 
the constraint condition rj\ X , N >= we find that the condition -^rj\x, N >= generates 
a new constraint operator which we identify as the eigenvalue constraint operator E given 
by[v,H] = £E: 

; E\ X ,N>=0. (3) 

x=0,L 

The Current Constraint 

The current constraint is obtained from the Noether current. We perform an infinitesimal 
gauge transformation C^(x)\0 >= e K ^C\x). Due to the folding this transformation trans- 
lates into: Cf(a;)|0 >= e ie ^C\(x)\0 > and Cj(a;)|0 >= e ie( --^Cl(x)\0 >. The state \x,N > 
must be invariant under a periodic gauge transformation e(x) = e(x + L). As a result of the 
transformation the Hamiltonian h = ^[^a,i(~ id x — ^-^i) 2 + 5 ai2 ( — idx + ^-fy) 2 ] is replaced 
byh= Jf-[6 aA (-id x -f <p 1 + d x (e(x))) 2 + 5 a , 2 (-id x + ^ ^ 2 + 9_ :r (e(-x))) 2 ]. The constraint is 
invariant under the gauge transformation vi\x)'q{x) = rf ! (x)rj(x). The constraint operator n 
is replaced by the transformed one rj = [e~ ie ^r](x)]\ x=0tL = \e~ lelyX 'Ci(x)— e^^C^fx)] | x=0 ,l 
, rj\x, N >= 0. [e(x) is an arbitrary periodic function in L, which is continuous at x = 
and has a continuous derivative d x (e(x)) ^ at x = 0. For example, any function with 
the Fourier expansion e(x) = Xlr=r ^rsin[^-x] and Fourier components ^7^0 obeys 

this condition.] The transformed constraint rj\x, N >= must hold at any time, therefore 
we have the equation : -^rj\x, N >= 0. Applying the Heisenberg equation of motion for 
the transformed Hamiltonian h and keeping only first order terms in d x (e(x)) that obey 



E 



2tt 



id x - — ipA d(x) 



2tt 

-idx + -j^^P2 ) C 2 (x) 
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d x (e(x))\ x=0 ^ gives us: 

.2 rL 



th-v\x, N >=^J o dx fc C\(x)(-id x - —ft + d x (e(xWCi(x) 
+Cl(x)(-td x + ^-0 2 + d_M-x))) 2 C 2 {x)]\ X , N >= 0. (4) 
Using the energy constraint E\ X , N >= we identify the current continuity constraint (3: 



P 



-id x - -j^r ) Ci(x) + ( /<),. + -j^02 } C-A/) 



\ X=0>L ; (3\ X ,N>=0. (5) 



To conclude the eigenstate \x, N > for N particles in two rings must satisfy the following 
equations: 

H\x,N>=E(N)\ X ,N> ; n\ X , N >= ; E\ X , N >= ; P\ X ,N>=0 (6) 

The eigenfunctions will be given in terms of the amplitude wave functions: For example 
the single particle state is given by : 

\ X ,N=1 >= J Q L dx[h{x)C\{x) + f 2 (x)Cl(x)]\0 > . 
Similarly the two particle state is given by: 

\ X ,N = 2 >= [ L dx [ L dy\f 1A (x,y)C!l(x)Cl(y)+f lfl (x,y)^ >]. 
Jo Jo 

(7) 

The amplitudes f\{x)J 2 {x) and fi t i(x,y), fi >2 (x,y) f2,2(x,y) are determined by the con- 
dition given in Eq. (8). 
The Current Operator 

The iV particle wave function < xn, ...Xi\x,N > must obey periodic boundary conditions: 
< 0\C ai (x 1 )..C olk (x k )..C aN (x N )\x, N >=< 0|C ai (xi)..C afc (a; fe + L)..C aN (x N )\ X , N > 

where ojj takes two values ccj = 1 or cti = 2. 

Once the eigenfunction \ X ,N > is known we can use the current operators J±(x) and 

J 2 (x) in the second quantized form to compute the current in each ring: 

Mx) = -^-[ct(x)(d x - i^<h)c x {x) - m - t^^ciix))^) ■ j x { X ) = <N : x y i \ x)l * ,I * > , 

%lm L L <JM, X \Xi^ > 

ti*) = ^\Pl(x)(d x + i^fr)C 2 (x) - ((d x + i ^ 2 )Cl(x))C 2 (x) ; J 2 (x) = < ^y^^t8 .) 
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II. THE EMERGING CONSTRAINT CONDITIONS FROM THE TIGHT BIND- 
ING FORMULATION 

The Hamiltonian in Eq. (1) must be supplemented by the coupling term between the 
rings. The most general form for the coupling is given by: 

H couplmg = [ dx6(x)[-U ± (& 1 (x)C 2 (x) + Ci(x)C 1 (x)) + U\\(Cl(x)C 1 (x) + 
Jo 

L 

dxS(x)[U ± (C\(x) - Cl(x))(C\(x) - C 2 {x)) + {U\\ - U ± )(Cl(x)d{x) 
+Cl(x)C 2 (x))}. (9) 

We introduce the notation U± = tU and U\\ = sll where s and t are parameters. Using the 
spinor representation C(x) = [C\(x) , C 2 (x)] T we can rewrite the coupling Hamiltonian in 
terms of the Pauli matrix o~\ and the identity matrix /: 

Hcoupiing = J L dx5(x)U^ + (x)[(sl - ta 1 )]^(x). 
This problem belongs to the class of delta function potential considered in Quantum Me- 
chanics. 

A-The wave function for a single particle, N=l 

\ X ,N=1 >= J Q L dx[h{x)C\{x) + f 2 (x)Cl(x)]\0 >, 

(H + H couplmg )\ X ,N=l >=E(1)\ X ,N=1 >. 
As a result we obtain the Schrodinger equation in terms of the two amplitudes fi(x), f 2 (x). 
The Hamiltonian in Eq. (1) together with H coup u ng can be solved using the method for delta 
function potentials. We integrate the single particle Shrodinger equation around x = 0, L 
and obtain the discontinuity derivative of the spinor ty(x) = [fi(x), f 2 (x)] T which obeys 

dV(x) i x=e _ dV(x) \ x =e 
dx \x=—€ — dx \x=L-e- 

[(-id x - ^)h{x)\lzl_ e = =§*U\[{sh{e) - tf 2 (e)) + (sf^L - e) - tf 2 (L - e))], 

(~id x + 4 f V2)f2{x)\lzl_ e = =§±U\[{sf 2 {e) - th{e)) + (sf 2 (L - e) - tj\(L - e))]. 

This set of equations gives us the boundary conditions for the present problem. Indeed 
these equations are determined by the discontinuity function U[sf 2 (0) — £/i(0)]. For this case 
the solution follows from the method of the delta function potential-see Griffiths Quantum 
Mechanics section, 2.5 page 73. 
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B-The wave function for two particles, N=2 

In order to compute the wave function for N particles we have to com- 
pute the boundary conditions for the amplitudes of the wave function. We will 
consider the case of two particles which can be generalized to many particles. 
\ X ,N = 2 >= j L d Xl j L dx 2 [f 1A (x 1 ,x 2 )Cl(x 1 )C\(x 2 ) + f 1 , 2 (x 1 ,x 2 )Cl(x 1 )Cl(x 2 ) + 
f 2 , 2 ( Xl ,x 2 )Cl( Xl )Cl(x 2 )}\0 >}. 

Using the eigenvalue equation: (H + H coupHng )\x, N = 2 >= E(2)\x,N = 2 > we inte- 
grate the two particle Shrodinger equation around X\ = 0, L and obtain the discontinuity 
derivative for the three amplitudes x 2 ), f 2 , 2 (xi,x 2 ), fi, 2 (xi,x 2 ) 

[(-id xi - ^)h,i{x u x 2 )\7=l_ t = ^UlMfafa = e,x 2 ) + f 1 , 1 (x 1 = L-e,x 2 ))], 
[(-id Xl - f ^i)/ 2)2 (a:i, x 2 )\ x x \Zl_ e = =§-Ul[-t(f 2 , 2 ( Xl = e, x 2 ) + f 2>2 ( Xl = L - e, x 2 ))}, 

[(-id xl -%$i)fiAx,X2)\Z\=L- e = ^ 
Similar equations are obtained by exchanging x\ with x 2 . This set of equations determines 

the two particle wave function < xi,x 2 \x,N = 2 >. This procedure is rather involved but 

can be generalized to the iV particles case. 

C-The strong coupling limit U — > oo 

Next we investigate the strong coupling limit and show that the problem can be simplified 

to a constraint problem. We consider the case s — t — 1 for which we have the scattering 

matrix S given by: 

S = Texp-^n^ <ir[(c 1 +(,=o,.)-c+(,=o,.))(c 1 (,=o,.)-c 2 ( a; =o > .))]_ FoT u ^ ^ the sca t ter ing 
matrix S obeys : 

Z; m ^re-^nodT[(c+(x=o,r)^^ >^ (Ci(a . = Q t) _ 

C 2 (x = 0,t))\ X ,N>=0. 

As a result the field (Ci(x = 0, r) — C 2 [x = 0, r)) is enforced to satisfy C±(x = 0, r) — C 2 [x = 
0, r) =0 which is equivalent to the constraint condition: 
T]\X, N >= [Ci(x) - C 2 (x)]U =0 ,l|x, ^ >= 0. 
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III. COMPUTATION OF THE WAVE FUNCTION FOR EQUAL FLUXES 



For the strong coupling limit we will use the constraints given by the equation 6. When the 
fluxes are the same for both rings the constraint operator /3 is simplified to a new constraint 

7 = = <p 2 ): 

7=[d x C 1 (x)+d x C 2 (x)]\ x=0iL ; 7 | x ,iV>=0. (10) 
The N particles wave function for equal fluxes must satisfy the following conditions : 

H\x,N>=E(N)\ X ,N> ; V \ X , N >= ; E\ X , N >= ; 1 \ X ,N>=0. (11) 

A-The single particle case 

The single particle case corresponds to one electron in two rings. The state for one 
particle is given by: \x, N = 1 >= J Q L dx[fi(x)Cl(x) + f 2 (x)C 2 (x)]\0 >. The two-component 
spinor amplitudes fi(x) and f 2 (x) represent the wave function. Using the Hamiltonian given 
in Eq. (1) we can write down the eigenvalue equation H\x, N = 1 >= E(1)\ X , N = 1 >. A 
standard calculation shows this equation is equivalent to two eigenvalue equations for the 
amplitudes fi(x) and f2(x). 

^(-*d x -^J f 1 (x)=E(l)f 1 (x) ; ^(- idx -^y f 2 (x)=E(l)f 2 (x). (12) 

The constraint operators given in Eq. (8) generate the followings boundary conditions at 
x = 0: 

/i(x = 0) = / 2 (a; = 0); [d x h{x) + d x f 2 (x)]\ x=0 = 0. (13) 

The first equation is equivalent to the continuity of the wave function at x = and the 
second equation describes the continuity of the derivative of the wave function (once we 
fold back the space ) at x — 0. From the eigenvalue equation given in Eq. (9) we find: 
E(n;N = 1) = £(f )\n - $f for the ring one and E(-n;N = 1) = £(f ) 2 (n - (p) 2 
for the second ring. Due to the folding of the space around x=0, the eigenvalue with the 
quantum number n in ring one and Quantum number — n in the second ring are equal. 

This result holds for the quantum numbers, n = 0,±1,±2, The single particle state 

\n, N = 1 > for (p 7^ \ is given by: 

\n- N = 1 >= -L f dx[e i2 -r nx C\{x) + e^™Cl{x)\\Q > . (14) 
v2L Jo 
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To understand this result we fold back the ring such that x — > —x. This means that if the 
particle in the first ring (x < 0) has the momentum it will be perfectly transmitted to 
the second ring with the same momentum and the same amplitude. If we remove the point 
x = and create a ring of a double length 2L, the current will be the same as in one ring 
with the same flux. Indeed, the only difference being the doubling of the size. As a result, 
we will have half of the current in a single ring. (If we rescale the length, we find the same 
current as in one ring [11].) It is important to remark that the states \n; N = 1 > and 
| — n; N = 1 > correspond to two different eigenvalues. Therefore, for a given eigenvalue we 
can not have a linear combination of waves e %2 f nx and e ~ l? r nx in the same ring. The wave 
e l — nx in ring one will be transmitted into the second ring without any reflection, the form 
of the transmitted wave will be ^r % — nx (hi the unfolded coordinates the form of the wave 
will be e l ~ ny in the second ring for y < 0). In Fig. 1 we show the current flow for two rings 
with equal fluxes in the unfolded geometry. The current vanishes if we have the opposite 
flux in the two rings, as depicted in Fig. 2. 

The case (p— \ deserves special consideration. The eigenvalue operator E has two pairs 
of momentum with the same eigenvalue: The first pair ri\ = n in the first ring and n 2 = —n 
for the second ring and the second pair n[ = —n + 2(p (ring one) and n' 2 = n — Kp (ring 
two). As a result we obtain two degenerate eigenstates \n; N — 1, + > and \n; N = 1, — > 
given by: 

\n; N=l,+ >= -L= I dx[e l2 -r nx C\{x) + e^™C\{x)\\Q > ; 
v2L Jo 

|n; N = 1, - >= -L [ L dx[e- i2 -r {n "^ )x C\{x) + e^ n ~^ C%{x)\\Q > . (15) 
v2L Jo 

As a result the current for the state \n; N = 1, — > will be opposite to the current for 
the state \n; N = 1,+ > Since the two eigenstates \n; N = 1,+ > and \n; N = 1,— > 
are degenerate, the single particle state will be given by two linear combinations of the 
eigenstates \n; N — 1, + > and \n; N — 1, — >: \x(n), (p — |; N — 1 >= a + \n; N — 1, + > 
±a^\n; N — 1, + > with the condition \a + \ 2 + |ct-| 2 = 1. For the special values |a + | 2 =|o;_| 2 
the current will vanish. 

B-The Two particles state 

We will construct the two particles state and will show that due to the constraints not 
all the antisymmetric combination of the single particles states which obey the constraints 
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are allowed. Imposing the constraints on the two particles state imposes further restrictions. 
The two particles eigenstate is determined by the three components fu(xi, X2), /12OE1, £2) 
and f 22(^X1, x 2 ) that obey the eigenvalue equations: 

fi[H^ - 2 f 0) 2 + (~id X2 - ¥0) 2 ]/n(xi,s 2 ) = E(2)f 11 (x 1 ,x 2 ), 
- 2 f 0) 2 + (id X2 - 2 fmfu(x 1 ,x 2 ) = E(2)f 12 (x 1 ,x 2 ), 

^[(id Xl - 2 f0f + (id X2 - 2 f£) 2 }f22(x 1 ,X 2 ) = E(2)f22(x 1 ,X 2 ). 

The amplitudes fu(xi, x 2 ), fi2(xi,x 2 ) and f22i.X1.X2) are constructed from the sin- 
gle particle states which are represented in terms of the complex coordinate Z[x) = 
e l2 ^ x and Z*(x) = e~ l2 ^ x . We introduce the antisymmetry operator A, which acts 
both on the space coordinates and the ring index matrices An (two particles on ring 
one), Au (one particle on ring one and the second on ring two), and A22 (two par- 
ticles on ring two). When the operator A acts on a two particle wave function 
it gives : A[A 1 2(Z(x 1 )) n (Z(x 2 )) m } = [A 12 Z(x 1 )) m (Z(x 2 )) n - A 21 Z(x 2 )) n (Z( Xl )) m ] and 
A[A u (Z( Xl )) n (Z(x2)) m } = [A ll Z(x 1 )) m (Z(x 2 )r - A u Z(x2)) n (Z(x2)) m } for i = 1,2. 

From the eigenvalue constraint E\n,m;N = 2 >= we obtain the condition for the 
eigenvalues. The only possible solution for these equations are states with m = —n which 
give eigenvalues E(2) = E(n, -n; N = 2) = ^(^) 2 [{n - fi) 2 + (-n - 0)\ n = 0, ±1, ±2... 

For amplitude fu(xi, x 2 ) we consider only the single particle states with n and —n which 
have the eigenvalue E{2) = ^(t") 2 ^ — + ( — n ~ l f > ) 2 ]- ^ e construct the antisymmetric 
amplitudes will be given by: 

fn(x 1 ,x 2 ) = A 11 [(Z(x 1 )) n (Z(x 2 ))- n - (Z(x 2 )) n (Z(x 1 ))- n }. 
Similarly for two electrons on the second ring f 22(^X1, x 2 ) we have: 

f22(x 1 ,x 2 ) = B 11 [(Z(x 1 ))- n (Z(x 2 )) n - (Z(x 2 ))- n (Z(x 1 )) n }. 
The amplitude for one electron on ring one and the second electron on ring two is given by 
fi2{x\, x 2 ) .This corresponds to two pairs of states n, n and — n , — n. The eigenvalue for the 
pair n , n is equal to E{2) = ^(^ ) 2 [{n — (p) 2 + (—n — (p) 2 }. For — n , — n we have the same 
eigenvalue. The amplitude /i 2 (xi,z 2 ) is given by the linear combination of the two pairs . 
Using the antisymmetry operator A we obtain the amplitude /i 2 (xi,x 2 ) for the two pairs: 
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f 12 ( Xl ,x 2 ) = A[A 12 (Z(x 1 )) n (Z(x 2 )) n }+A[B 12 (Z(x 1 ))- n (Z(x 2 )r n ] 
= [A 12 (Z( Xl )) n (Z(x 2 )) n - A 21 (Z(x 2 )r(Z( Xl )) n } + 
[B l2 {Z{ Xl ))- n {Z{x 2 ))- n - B 21 (Z(x 2 ))- n (Z( Xl ))- n ]. 

Using constraints given in Eq. (8) for the two particles state \n, m; N = 2 >: rj\n, m; N = 
2 >= 0, E\n,m;N = 2 >= and j\n,m;N = 2 >= 0, and we obtain the following 
boundary conditions: 

2/n(a;i,0) = /i,2(zi, 0) ; [2d x J n (x u x 2 ) + d x J 12 (x u x 2 )] X2=0 = 0, 

2/2,2(2:1, 0) = /i, 2 (0, xi) ; [2d X2 f 22 (x 2 , xi) + d X2 f l2 (x 1 , x 2 )] X2=0 = 0. 

From these equations we find that the amplitudes obey the relations: Ai 2 = —A 2 \ = 2 An; 
An — A 22 and B 2 i = —B 12 = 2A 22 . We introduce the antisymmetric spinor notation 
6i ; 2 = ^r, which obeys the relations: e\l^ = — e 2-i an d { e \^Y ' e i|2 = 1 (^ ne upper index 
1, 1 means that we have one electron in each ring, the bottom index 1,2 or 2, 1 represents 
the order . e}|2 the first electron is ring one and the second electron is on ring two and e^.'i 
represents the first electron on ring two and second electron on ring one.) The normalized 
two particle state is given by: 

|n 5 -n;N = 2 >= 

[ dxi [ dx 2 [fii(x u x 2 )C\(xi)Cl(x 2 ) + fi 2 {xi,x 2 )C\(xi)Cl(x 2 ) + f 22 (x u x 2 )Cl(xi)Cl(x 2 )\\Q > 
Jo Jo 

= [ dxi [ dx 2 ^-[[(Z( a : 1 ))™(Z*(x 2 ))™-(Z( a : 2 ))"(^(x 1 ))"]C 1 t ( a ; 1 )C 1 t ( a ; 2 ) 
Jo Jo 4i/ 

+2e 1 ^[Z(x 1 )r(Z(x 2 )r - {Z*{x 2 )nZ*{xi)T]C\{xi)Cl{x 2 ) 

+[(z*( Xl )nz(x 2 )r - (z^x 2 )nz(xi)r]ci(xi)ci(x 2 )]\o >\. ae) 

The off-diagonal spinor component fi 2 (xi, x 2 ) oc 4ism(^n((xi + x 2 )) is symmetric in space 
and resembles the BCS pairing wave function (once we identify the ring index with the spin) 
in contrast to the diagonal elements fn((xi, x 2 ) and f 22 {{xi, x 2 ), which are antisymmetric 
in space. The two particles state, which obeys the constraints are different from the two 
particles state constructed from the single particles, which obey the constraints. Using the 
single particle states \n; N = 1 > and \m; N = 1 > [which obey Eq. (11)] we construct an 
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antisymmetric tensor product \n, m; N = 2 >buiid— \n; N = 1 > \m; N = 1 > —\m; N = 1 > 
\n; N — 1 >. This state is not a solution which obeys the constraints for the two particles 
state. The only possibility is to have an antisymmetric tensor product of two states with 
vanishing total momentum \n, —n; N = 2 >= \n; N = 1 > | — n; N = 1 > — \ — n; N = 1 > 
\n; N — 1 >. (The ground state for the two particles {(p < |) is given by the eigenstate 
|1, — 1; N = 2 > .) This structure persists for an even numbers of electrons N = 2M and 
gives rise to a robust state absent for the single ring. 

C-The three particles state 

The wave-function for three particles can only be found for special configurations 
| to, n, — n; N = 3 > m ^ n and m ^ —n. The ground state will be given by the 
state |0, 1,— 1;N = 3 >. The three particles state is determined by the four amplitudes 
fui(x 1 ,x 2 ,x 3 ), fu 2 {x 1 ,x 2 ,x 3 ), fi22(xi,x 2 ,x 3 ) and f 122 (xi , x 2 , x 3 ) , which obey the eigenvalue 
equation: 

^[(-id Xl -^f + {-id X2 -^fH-id Xz ^ 

£ i [(-id xl -^) 2 +(-id X2 -^) 2 +(-id X2 + 2 £^ 
^\{-id xl -^) 2 +(-id X2 +^y+{-id x .^^ 

£[Hd x + 2 f& 2 + (-id X2 + 2 f0) 2 + (-id X3 + f 0) 2 \f^{x u x 2 , x 3 ) = E(3)f 222 (x 1 ,x 2 , x 3 ). 

Using Eq. (8) we obtain the following relations for the spinor components: 

3/ m (a;i,X2,0) = fn 2 (x 1 ,x 2 ,0) ; [3d X3 f lu (x 3 , x u x 2 ) + d x J 112 (x 2 , x u x 3 )] X3=0 = 0, 
3/222(0, xi,x 2 ) = /i 22 (0,xi,x 2 ) ; [3d x3 f 222 (x 3 ,xi,x 2 ) + d z f U2 (x 2 ,x 1 ,z 3 )] Z3= o = 0, 
2f 121 (x 1 ,x 2 ,0) = fi 22 (x 2 ,x 1 ,0) ; [3d X3 f 121 (x 1 ,x 2 ,x 3 ) + d X3 f 122 (x 1 ,x 2 ,x 3 )] z=0 = 0. 
The solution of the constraint equations fixes the eigenvalue and the state. The ground 

state eigenvalue is given by E g (0, 1, -1; N = 3) = ^(^Y[(0Y + (1 - 0Y + (-1 - vY) and 

the three particles ground state is : 

|0,1,-1;JV = 3 >= 

i*L pL pL 

/ dxi \ dx 2 \ dx 3 [$ ,i,-i(x 1 ,x 2 ,x 3 )Cl(x 2 Cl(x 2 )Cl(x 3 ) 
Jo Jo Jo 

+ 3 [ e U;2 ^2 Pi,x 3 (®0,l( x U X 2 )Z(x 3 ) - $ ,-l(^l, X 2 )Z*(x 3 )) + $ ,1 ,-l(^l, X 2 , X 3 )]C\{xi)C\(x 2 )Cl(x 3 ) 

+3[e}; 2 , 2 (^0,1(^2, x 3 )Z( Xl ) - $0,-1(^2, x 3 )Z*( Xl )) + $0,1 -ifa, x 2 , x 3 )]C}(a;i)Cj(j/)Cj(a;3) 

i=xi,X2,x% 

+% :1 ,_ 1 (x 1 ,x 2 ,x 3 )Cl(x 1 )Cl(x 2 )Cl(x 3 )]\0 > . (17) 
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This state is expressed in terms of the Slater determinants for two and three particles 
&o,±i(xi,x 2 ) , $o,i,-i(xi,x 2 ,x 3 ). (Pi,x 3 is the coordinates interchange operator defined by: 
P i)Z F(x l , x 2 ; x 3 )F(x 1 , x 2 ; x 3 ) = 5 i)Z F(x u x 2 ; x 3 ) + 8 ijX1 F(x 3 , y 2 ; £1) + 5 i:X2 F(x 1 ,x 3 ; x 2 )). The 
three particles states can be rewritten as an antisymmetric tensor product of the three single 
particles states, which obey Eq. (11): 

|0, 1, -1; N = 3 >= Ep(-1) P |0p(i); N=l> |1 p(2) ; N=l>\- l P{3) ; N = 1 >. 

D-The four particles state 

The wave function for four particles has the structure \n, —n, m, —m; N = 4 > 
with n 7^ m. The ground state is given by: 1 1, — 1, 2, — 2; AT = 4 > with the eigen- 
value E g (l, —1, 2, —2; N = 4). From Eq. (8) we find: H\l, —1, 2, —2; N = 4 >= 
£(4)|1,-1,2,-2;7V = 4 >, 77 j 1 , —1, 2, —2; N = 4 >= 0, E\l, —1, 2, —2; N = 4 >= 
and 7|1, — 1, 2, — 2; TV = 4 >= we obtain a set of equations for the spinor com- 
ponents fuu(x 1 ,x 2 ,x 3 ,x 4 ), fiu 2 (x 1 ,x 2 ,x 3 ,x 4 ),f 1122 (x 1 ,x 2 ,x 3 ,x 4 ) Ji 222 (x 1 ,x 2 ,x 3 ,x 4 ) and 
f 2222 (x 1 ,x 2 ,x 3 ,x 4 ): 

4:fuu(xi,x 2 ,x 3 ,0) = /ni2(^i, ^2,^3,0) ; 

[49 I4 / m i(xi,X2,X3,a;4) + ^ 4 / m2 (xi, x 2 , x 3 , x 4 ))} X4=0 = 0, 

4/2222(^1,^2,^3,0) = -f 1222(0, x 1 ,x 2 ,x 3 ); 

[4d x J 2222 (x 1 ,x 2 ,x 3 ,x 4 ) - d X4 f 1222 (x 1 ,x 2 ,x 3 ,x 4 )] X4=0 = 0, 

3/1112(^1, x 2 , 0,x 4 ) = -2f 1122 (x 1 ,x 2 ,x 3 ,0); 

[3^/1112 (x 1 ,x 2 ,x 3 ,x 4 ) -2d x J 1122 (x 1 ,x 2 ,x 3 ,x 4 )] X4=Q = 0, 

3/1222(^1,^2,2:3,0) = -2/1221(^1, x 2 , x 3 ,0); 

[^d x J 1222 (x 1 ,x 2 ,x 3 ,x 4 ) + 2d x J 1221 (x 1 ,x 2 ,x 3 ,x 4 )] Xi=Q = 0. 
The eigenvalue and the eigenfunction are: 
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E g (l, -1, 2, -2; iV = 4) = ^(^) 2 [(1 - y3) 2 + (-1 - ^) 2 + (2 - 0f + (-2 - <p)\ 
|1,-1,2, -2;iV = 4 >= 

/»L /»L /*L /»L 

/ dxi / drr 2 / rfx 3 / rfx 4 [$i-i, 2 -2(^i,^2,a:3,X4)C' 1 t (a;i)C' 1 t (x2)C'i t (a;3)C' 1 t (x4) 
Jo Jo Jo Jo 

+ 4e lil;2[ Pi,x 4 [®2,l,-l(xi,X 2 ,X 3 )(Z(x 4 )) 2 ~ § , X 2 , X 3 )(Z* (x 3 )) 2 

i=Xl,X2,X3,X4 

+$1,2,-2(31, x 2 , x 3 )Z(x 4 ) - $_ lj2 - 2 {x u x 2 , x 3 )Z*(x i )]]C\{x 1 )C\(x 2 )C\(x 3 )Cl(w) 
+6e?;i ;2 , 2 [[ J] X A,x 4 ] [$1,-1(31,32)^2,-2(33, 34) 

i=Xl,X2,X3 i=Xl,X2,X4 

+$1,2(31, x 2 )$_! _ 2 (x 3 , x 4 )]]C'I(xi)C}(x2)Cj(a;3)C , 2( a; 4) 

+ 4e !;2,2,2[ X [$2,1,-1(32, x 3 ,3 4 )(Z(:ei)) 2 - $_ 2 ,i,-i(x 2 , £3, £ 4 )(Z*(xi)) 2 

+$1,2,-2(32, 3 3 , x i )Z(x 1 ) - $-1,2,-2(32, 3 3 , x 4 )Z*(xi)]]C 1 t (a; 1 )C' 2 t (x2)C 2 t (a;3)C 2 t (a;4) 
+$1,-1,2,-2(31, 3 2 ,x 3 ,x 4 )C 2 t (a;i)C 2 t (a;2)C2 t (33)C2 t (34)]|0 > 

= X(-1) P |1 P(1) ; N = 1 > I - 1 P(2 ); AT = 1 > \2 P{3) ;N = 1 > | - 2 P(4) ; JV = 1 > . (18) 
p 

Where $1-1,2,-2(^1, 3 2 , x 3 , x A ) , $±2,1,-1(^1, x 2 , x 3 ) and $„, m (xi,x 2 ) are the S7ater 
determinant for 2, 3 and 4 particles. Here ef ^ x 2 an d e i^i;2,2 are the antisymmetric ten- 
sors for the ring index. 

E- The 2M particles state 

The 2M particles state is built from the single particles states n±, ..nk, -.nu given by Eq. 
(11) with vanishing total momentum: 

|ni, -n 2 , ..n 2 fe_i, -n 2k , ..n 2M -i, -n 2M ; N = 2M >= 

X(-l) p |n P(1) ; N =1>\- n P{2) ; N = 1 > ..> p(2 m-i); N = 1 > | - n P(2M ); AT = 1 > . 

(19) 

The ground state and the ground state energy are: |1, — 1, ...M, — M; AT = 2M > s = 
E P (-l) P |lp ( i);iV = 1 > I - l P(2) ;Ar = 1 > |2 P(3) ;Ar = 1 > I -2 P(4) ;AT = 1 > 
••••|fcp(2fc-i) ; jv=i > I - k P(2k)] N = 1 > ...|M P(2M -i); N =1>\- M P{2M) ; N = 1 > ; 
E a (l, -1, .., k, -k, ...M, — M) = * f Ef =1 P - V? + (-* - ^) 2 ]- 

F- The Current for equal fluxes 
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The current for equal fluxes with 1, 2, 3, 4 and 2M particles is the same in both rings: 

j N =i _ <N = 1; n\Ji(x)\n; N = 1 > _ H 2ir Ap-n _ 
Jl ~ <N = l;n\n;N=l> ~ [ m L l[ 2L 1 ' " " ±ij ±2 " 

T N=i ( , lx <N=l;<p=±, X {n)\j 1 {x)\x{n),<p = ±-,N = l > h 2n 2 2 
Jl ( ^ 2> = <^ = l;# = I >X (n)|x(n)^ = i^=l> = [ mT^" " I "" ][ 



tAT=2 _ 



< N = 2 



< AT = 2; —1, 11, — l;iV = 2 > m L 2L 



^=4 _ < iV = 4 



-1, 1| Ji(x)|l, — 1; TV = 2 > _ h 2tt 2<£ 

L™ r J Lr> r J 



-l,l,0|Ji(x)|0, ,1,-1; JV = 3 > _ ft 2?r 3<£ 

L™ r J Lr> r J 



<iV = 3;-l,l,0 0,,l,-l;JV = 3> l m L sl 2L s 



-2,2,-l,l|Ji(a;)|l,-l,2,-2;JV = 4 > _ 2tt 4y3 

L™ r J Lr> r J 



< JV = 4;-2,2,-l,l|l,-l,2,-2;JV = 4 > L m L"2L J 
jv =2M _ < iV = 2M; — M, M, ... - 1, 1| Ji(x)\l, -1, ...M, — M; AT = 2M > g _ h 2ir 2M(p 
1 " <N = 2M; -M, M, ... — 1, 1 1 1, -1, ...M, — M; N = 2M > g ~ ml" 2L 

The magnetization is given by the current area product: = 2Jf^-. For an 

even number of electrons we find that the current in a single ring is twice the current in a 
double ring J^ngie^-Hng = 2Jf =2M . The factor of \ is a result of the two component spinor 
state renormalization. At finite temperatures the two rings excited states have the form : 
|1, —1, ...M + p, — (M + p); N = 2M > e where p are integers. This state carries the same 
current as the ground state |1, — 1, ...M, — M; N = 2M > g . Therefore, we conclude that 
for an even (fixed) number of electrons the current will be the same at any temperature! 
(When the total number of electrons fluctuates, N N ± 2 thermal effects will decrease 
the current.) The situation for the odd number of electrons is different. Even for the two 
states |1, -1, ...M, -M, n = (M + p); N = 2M + 1 > and |1, -1, ...M, -M, n = — (M + 
p); N = 2M + 1 > we have different eigenvalues and at finite temperatures these states 
carry a different current. Therefore, the total current carried by all the states will be 
reduced like we have for a single ring where the unrestricted structure of the wave function 
allows any configuration of momenta, which generate an antisymmetric wave function in 
space: f i - sm9le ' rm9) {xi,x 2 , ...x n=2 m) = $n 1 ,n 2 ,....n 2M {xi, x 2 , ...x n=2 m)- To probe this even-odd 
structure we propose to attach a gate voltage to the rings. As a result, the magnetization 
will vary with the varying gate voltage. 
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IV. THE WAVE FUNCTION FOR OPPOSITE FLUXES 



For this case the constraint operator 7 is modified to: 7 = [{—id x — ^(p)(Ci(x) + 
C 2 {x))]\ x= o t L. For the single particle case we find the following boundary conditions: 

fi(x = 0) = f 2 (x = 0) ; -i[d x h(x) + d x f 2 (x)]\ x=0 ,L = f0\fi(x) + f 2 (x)]\ x=0 , L . 
We find that for this case the wave function must vanish. Only for integer values of flux 
n = integer = (p we have finite solutions fi(x) = f 2 (x) = e l ^ x with a vanishing persistent 
current. This result is in agreement with the fact that at the common point between the 
rings the fluxes must satisfy (p 2 = <fi + n. Therefore, the boundary condition can be satisfied 
for this case only if the the wave function vanishes at the common point . We mention that 
for two separated rings threated by opposite fluxes the magnetization will be zero only at 
the symmetry points. This result allows to control the current in one ring by reversing the 
flux in the second ring. 



V. TWO COUPLED CYLINDERS 



In order to build a theory which can be compared with the experiment we have to consider 
effects of interactions and effects of finite width geometry [13]. For realistic considerations 
the point contact between two rings is replaced by two narrow cylinders of height d <^ L 
which are in contact at the point (x = 0,0 < z < d The gluing condition is implemented by 
two narrow cylinders of height replacing the constraints in Eq. (6) by n(z, x = 0)\x, N >= 
and 7(2;, x = 0)\x, N >= 0. In the absence of disorder we obtain for each transversal channel 
r = 1, 2..r max one dimensional constraints: r) r (x = 0)\x, N >= 0; 7 r (x = 0)|x, N >= 0. 

Therefore, the current in the channel r is the same as the result given in Eq. (17). For 
iV electrons the current will be determined by the partition of iV electrons in the different 
channels : N = N x + N 2 + ..N r ... + N Tmax 

In the absence of disorder the current in cylinder one, at T=0 will be given by: 

jN _ \h 27rir 2(A r i+A r 2+-A f r-.A r r ma3: )y i _ r h 27rir2jVyi 
J l —VmLlV 2L \ — Vm L\V 2L I' 



16 



VI. CONCLUSION 



To conclude, a new method for applying constraints has been presented. This method 
has been used to compute the wave function for coupled rings. For an even number of 
electrons, only states with total vanishing momentum are allowed giving rise to a large 
persistent current and magnetization. For odd number of electrons at finite temperature 
the current and the magnetization are suppressed. We propose to confirm this even-odd 
effect by attaching the two rings to a varying gate voltage. Reversing the flux in one ring 
will cause the current to vanish in both rings. We construct the many article ground state 
which obey the constraints and show that not all the many particle states which are build 
from single particle states which obey the constraints are allowed. 

1. K. Sasaki, Y. Kewazoe, and R. Saito, Physics Letters A 321, 369-375 (2005). 

2. H. Aoki, J. Phys. C. 18, 1885-1890 (1981). 

3. Paul A. M. Dirac, "Lectures on Quantum Mechanics," Belfer Graduate School of Science, 
Yeshiva University, New York, 1964. 

4. M. Henneaux and Claudio Teitelboim "Quantization of Gauge Systems", Princeton Uni- 
versity Press (1992). 

5. Y. Avishai and M. Luck, J. Phys. A: Math. Theory 42, 175301 (2009). 

6. S. Itoh and S. Ihara Phys. Rev. B. 48, 8323 (1993). 

7. H. Bluhm, N.C. Koshnick, J.A. Bert, M.E. Huber, and K.A. Moler. Cond-Mat/08104 

8. Y. Aharonov and D.Bohm, Phys. Rev. 115, 485 (1959). 

9. N. Byers and C.N. Yang, Phys. Rev. Lett. 7, 46 (1961). 

10. M.Buttiker, Y. Imry, and R. Landauer, Phys. Lett. A 96, 365 (1983). 

11. L.P Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, Phys. Rev. Lett. 64, 2074 (1990). 

12. D. Schmeltzer, J. Phys. Condens. Matter 20, 335205 (2008). 
13.M.Grochol et al. Phys.Rev.B 74, 115416 (2006) 



17 



